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Abstract 

We consider the wave and ray dynamics of the electromagnetic field in 
a parabolic dome microcavity. The structure of the fundamental s-wave 
involves a main lobe in which the electromagnetic field is confined around 
the focal point in an effective volume of the order of a cubic wavelength, 
while the modes with finite angular momentum have a structure that 
avoids the focal area and have correspondingly larger effective volume. 

The ray dynamics indicate that the fundamental s-wave is robust with 
respect to small geometrical deformations of the cavity, while the higher 
order modes are unstable giving rise to optical chaos. We discuss the 
incidence of these results on the modification of the spontaneous emission 
dynamics of an emitter placed in such a parabolic dome microcavity. 
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1 Introduction 

The miniaturization of optoelectronic devices such as light emitting diodes or 
semiconductor lasers, is expected to lead to an improvement of their energy 
efficiency and to a lowering of the lasing threshold. This tendency towards 
miniaturization has led to the exploration of optical microcavities whose dimen- 
sions are of the order of a few wavelengths [EJ . In such microcavities the extreme 
confinement of the electromagnetic field modifies the interaction of the active 
medium with the radiation field so that the process of spontaneous emission 
is altered both in its spatial and its dynamical characteristics. Spontaneous 
emission can thus be redirected, enhanced or inhibited in a way that may be 
exploited for the operation of light-emitting diodes or lasers. A modification 
of the characteristics of spontaneous emission, such as its directionality or the 
emission rate, has been shown for several microcavity designs such as for the 
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Figure 1: Atomic Force Microscope image of a "hill" of diameter 7.2 fim and 
parabolic cross-section of height 1.8 /an, etched on a GaAs substrate by a Fo- 
cused Ion Beam apparatus. When covered with gold it constitutes a concave 
parabolic mirror with its focal point inside the GaAs substrate. 

traditional Fabry-Perot planar cavities § and for disk-shaped || or spherical 
[f| cavities displaying whispering gallery modes. 

One of the key requirements for enhancing the dynamics of spontaneous 
emission and reducing the laser threshold is that the electromagnetic field at 
the site of the emitting dipole should be enhanced inside the cavity with re- 
spect to its value in free space. A class of resonators for which this can be 
achieved very efficiently are confocal cavities: A few experiments with spheri- 
cal confocal cavities ||, or semi-confocal microcavities || have been reported 
already, in which significant spontaneous emission modification or extremely 
low laser thresholds have been observed. Among the different designs of con- 
cave mirrors, parabolic mirrors have an important advantage in that their focal 
point displays no astigmatism and is free from spherical aberrations. Basic ge- 
ometric optic thus leads us to expect that double-parabolic confocal cavities or 
piano-parabolic semi-confocal cavities should display a strong enhancement of 
the electromagnetic field in the vicinity of the focal point, and a concomitant 
modification of the emission characteristics of an active medium placed there. 

This paper presents a theoretical analysis of microcavities formed by a 
parabolic mirror at or close to the confocal condition. The study is motivated 
by experimental work in which such a system has in fact been fabricated. The 
experimental characterization of the modal structure and dynamics, being now 
in progress, will be given in a separate publication [Q. Here, we briefly describe 
the experimental structure, in order to define the system for which our model 
calculations are intended. We have fabricated a semi-confocal piano-parabolic 
semiconductor microcavity (see Fig. hi) by etching an appropriately-prepared 
GaAs wafer by a Focused Ion Beam || to produce a "hill" of cylindrical sym- 
metry and parabolic vertical cross-section having a diameter of 7.2 /urn and 
a height of 1.8 fim (corresponding to optical lengths of respectively 27 A and 
6.75 A for a wavelength (in vacuo) of 960 nm) which was subsequently covered 
with a thin metallic layer of gold. This gold dome constituted thus a concave 
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Figure 2: Schematic cross-section of a semiconfocal parabolic dome cavity, con- 
sisting of parabolic gold mirror and a planar Bragg mirror placed at the focal 
plane of the parabola. The cavity spacer is made of GaAs, and the light emitter 
is a quantum well placed in the vicinity of the focal plane. 

parabolic mirror with its focal point inside the GaAs substrate. At the base of 
the parabolic hill, the wafer had a 6-period GaAs/AlAs Bragg mirror, closing 
the semi-confocal cavity (see Fig. |J). This cavity is expected to possess a mode 
in which the electric field is strongly enhanced in the vicinity of the focal point, 
so that a localized semiconductor emitter, such as a quantum box or quantum 
well, emitting at a wavelength near 960 nm, will have its spontaneous emission 
greatly enhanced when placed there. The use of a dielectric mirror with lower 
refractive index rather than a metallic mirror at the focal plane is important 
because it introduces a boundary condition that requires the tangential electric 
field to be maximal at the focal plane. This condition cannot be fulfilled on a 
metallic mirror, on which the tangential electric field should vanish, producing 
thus a vanishing field at the focal point of the parabola. 

In order to understand the operation of such a cavity and to assess its per- 
formance in modifying spontaneous emission, in this paper, we examine first the 
modal structure of an ideal confocal double-parabolic, or semi-confocal piano- 
parabolic microcavity. We then investigate the stability of these modes with 
respect to geometric deformations of the cavity that correspond to deviations 
from confocality; this condition is inevitably violated in a realistic cavity due 
to fabrication defects. The discussion of this case provides the conceptual and 
theoretical background for the experimental analysis to be presented in a sub- 
sequent paper. 

The calculation of the modal structure of the parabolic dome microresonators 
cannot be treated within the paraxial approximation of conventional |Ti"| ] res- 
onator theory, because of the very large aperture displayed by the parabola and 
because the cavity dimensions are comparable with the optical wavelength. Ex- 
tensions of the paraxial approximation to the highly convergent (or divergent) 
beams produced by parabolic mirrors are cumbersome even in macroscopic res- 
onators H where the optical axis is long compared with the wavelength - in 
microresonators, the latter breaks down as well. However, there are other ap- 
proximate techniques which are well-suited to the problem we consider. As 
a valuable tool for simplifying the exact solution of Maxwell's equations for 
the cavity modes, we employ a short-wavelength approximation leading to sim- 
ple WKB quantization conditions. The assumption that wavelengths are much 
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shorter than the relevant cavity dimensions is common to both WKB and parax- 
ial approximation, and it is therefore at first sight surprising that the WKB ap- 
proach yields excellent quantitative agreement with the exact cavity spectrum 
even for the longest-wavelength modes of the parabolic cavity. We show how 
this arises by discussing in detail the structure of the classical ray dynamics in 
the resonator which makes the WKB approximation possible. As a result, we 
shall then also be able to assess the stability of the modal structure with regard 
to fabrication imperfections, based on a ray analysis for parabolic cavities in 
cases where confocality is violated. To characterize the modes of the parabolic 
resonator, the internal caustic structure formed by the rays turns out to be of 
crucial importance. These considerations establish a connection between the the 
microcavity optics of the paraboloic dome and the field of quantum chaos: even 
minute deviations from confocality introduce chaos into the ray dynamics, and 
we have to address the significance of this effect for the relevant cavity modes. 

The paper is organized as follows: Section 2, introduces the mathematical 
model that describes confocal parabolic cavities, while Section 3 presents the 
wave equation for the electromagnetic field in cylindrical and parabolic coordi- 
nates and discusses its exact vectorial and scalar solutions. Section 4 presents 
the WKB approximation of the wave equations for the parabolic cavity, an ap- 
proach that will permit us in Section 6 to make the connection with ray optics, 
while Section 5 compares the numerical solutions of the wave equations in the 
parabolic microcavity with those of the WKB approximation. Section 6 intro- 
duces the main concepts of ray optics applied to our parabolic cavities with 
cylindrical symmetry, while Section 7 analyzes the stability of the ray trajecto- 
ries in a parabolic cavity in which there is a slight deviation from confocality. 
Section 8 discusses the problem of the finite acceptance angle of Bragg mirrors, 
a feature that limits the lifetime of modes in semiconfocal cavities bounded by 
such mirrors. Finally, Section 9 summarizes the results of this study and gives 
its conclusions. 



2 The model 

We consider a model structure for an ideal semi-confocal cavity which is bounded 
by a metallic concave parabolic mirror on one side and a planar dielectric mirror 
on the other side, placed at the focal plane of the parabola. 

In cylindrical coordinates (p, z, (f) the parabolic mirror is given by 

-i-V (" 

where / is the focal distance of the parabola, while the focal plane (and the 
planar mirror) corresponds to 

z = (2) 



5 



Figure 3: (a) Representation of parabolic cylinder coordinates in the z - p plane 
of a cylindrical coordinate system [z is the vertical axis). The third dimension 
is obtained by rotating the figure around z by the angle <f>. The point A is 
specified by £ = 1.3/, T] — 0.9/ and = 0. The focus of all parabolas is at 
the origin, (b) By unfolding the parabolic dome into a double paraboloid, the 
boundary conditions on the common focal plane can be restated as simple parity 
requirement under reflection at this plane (z = 0). For TE modes, the electric 
field must be even under this reflection. The unfolded cavity is shown in side 
view with meridians which make 90°-corners at the focal plane. The latter is 
also the equatorial plane of the cavity. 



It is convenient to describe this cavity in parabolic coordinates (£,77,0), whose 
properties are summarized in the Appendix. For reference, we reproduce here 
the transformation to cylindrical coordinates as given in Eq. J99]): 

{ P=VFl /OX 

To illustrate this coordinate system, we show in Fig. |^ (a) how the intersection 
of the coordinate surfaces defines a point A in the plane z versus p. Also shown 
is the cavity shape itself: the parabolic mirror corresponds to 

£ = 2/, (4) 

and the planar dielectric mirror is at 

i = n. (5) 

In an ideal cavity, the parabolic metallic mirror can be assumed to be lossless, 
displaying an amplitude reflectivity r = -1. This produces a n phase change upon 
reflection so that it corresponds to a boundary condition in which the tangential 
electric field vanishes. In parabolic coordinates this can be expressed as 

£„(£ = 2/) = E^ = 2f) = Q £? € (£ = 2/) = (6) 

Similarly, the planar dielectric mirror can be assumed to have a reflectivity of r — 
+ 1, producing no phase change upon reflection so that the tangential magnetic 
field vanishes on the focal plane of the cavity. In cylindrical coordinates, this 
can be expressed as 

B p (z = 0) = B^z = 0) = and E z (z = 0) = (7) 



G 



Alternatively, this implies that the tangential electric field is maximum on the 
focal plane and is symmetric under reflection of the whole cavity at the z = 
plane. Thus, instead of considering this plane as an additional boundary with 
the properties (^) , one can unfold the cavity across this plane by reflection, to 
obtain a confocal double paraboloid shown in Fig. || (b). This extended cavity 
requires only the metallic boundary conditions on its parabolic walls, that is 
Eq. (|^) and its equivalent in which £ and n are interchanged. It will support 
modes that can be either symmetric or antisymmetric under reflection at the 
focal plane. If we restrict ourselves to modes in which E p and are symmetric, 
this subset is identical to the modes of the original dome with the conditions of 
Eqs. © and (|). 

The advantage of considering the unfolded cavity is that the focal plane as a 
physical boundary drops out of the discussion; this will considerably simplify the 
interpretation in terms of the ray picture later on. Therefore, in the remainder 
of this paper, we can refer to Fig. [| (b) as our model system. 



3 Wave equation 

The electric field E obeys the vectorial wave equation 

d 2 E 

V x V x E + lie—- = (8) 
ot 

under the additional constraint that its divergence must vanish 

V--E = 0. (9) 

The boundary conditions and the constraint of zero divergence imposed on 
the electromagnetic field in general will lead to a coupling between the various 
vectorial components of the electric and magnetic fields. In simpler geometries 
such as cylinders, spheres or rectangular cavities, a suitable choice of polariza- 
tions reduces the problem to finding the eigensolutions of a scalar Helmholtz 
equation However, in our case the three polarizations and the intersecting 
parabolic surfaces forming the resonator cannot be labeled by the coordinate 
lines of a single orthogonal coordinate system, as is possible in the textbook 
systems mentioned. We discuss now the implications of this complication. 

3.1 Vector field components in cylindrical coordinates 

After combining Eqs. (||) and (Q) to the wave equation, 

V 2 E - /je^J = 0, (10) 
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we can take advantage of the cylindrical symmetry of the problem by expressing 
the wave equation for a time-harmonic electric field oscillating at frequency lo 
in cylindrical coordinates (p, z,(j)), as 

' V 2 E p -j lEp -j I ^+ neu 2 E p = (a) 
< V^-^ + ^ff + M ec^ = (b) (11) 

V 2 E Z + peu 2 E z = (c) 

We note that the wave equation couples the radial and angular components of 
the electric field (E p and E$), while the equation for the axial component E z is 
scalar. One can achieve a further simplification in this system of equations as 
follows: 

The rotational symmetry around the z axis permits us to assume a </>- 
dependence of all components of the field of the form 

Q(p,z)-e im t (12) 

With this ansatz, Eqs. ([ll] a,b) can be written as 

p 2 [V 2 + k 2 ] E p -E p = 2imE <t> (a) 

(13) 

p 2 [V 2 + k 2 ] E 4> -E <j} = -2imE p (b) 

where k = ^fjleuj is the wavenumber inside the parabolic dome. If the azimuthal 
quantum number m = 0, this reduces to two identical equations. If, on the other 
hand, m 7^ 0, we can form a suitable linear combination of E p and E^ which 
decouples these two equations. Naively setting E$ — would not achieve this 
goal because it forces both field components to vanish. 

The proper linear combination in which to decouple this system of differential 
equations is obtained with the definition 

E P = ^(E+-E_) 

(14) 

Et=j s (E + + E-). 

Then E± are the solutions of the equations 

p 2 [V 2 + k 2 ] E± = (l± 2m)E ± (15) 

This definition again makes use of the azimuthal symmetry of the resonator, 

which implies that the circular polarizations a± = T^tj (p± icj^j are decoupled 

in the cylindrical wave equation. In this way, we have therefore formally decou- 
pled the original system of equations Eq. (O) for the vector field components. 
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In the special case m — 0, case E+ and E— will moreover be linearly dependent 
because their respective equations again coincide. 

However, this decoupling of polarizations in the wave equation does not 
reduce the problem to a truly scalar one because the field components are still 
coupled by the boundary conditions and by the condition of zero divergence. 
On the "top" parabolic mirror, the conditions (^J) in terms of the cylindrical 
components of the electric field now read: 

( i^J{E + -E_) + ^E z =Q 



at £ = 2/ i E. 

drj 



(16) 



The first line expresses the condition E v = 0, the second and third lines represent 
E& = and = 0, respectively. On the "bottom" parabolic mirror the 
boundary conditions are the same as in Eq. ( |l6|) with £ and r\ interchanged. 



3.2 The absence of longitudinal electromagnetic modes 

Unfortunately, the set of boundary conditions Eq. ( |l6| ) is not yet a complete 
list of constraints that we have to satisfy. An additional requirement is that 
the field at every point in the resonator has to have zero divergence, which in 
parabolic coordinates reads 



E+ - E- + m(E+ +E_) + — .(— + —)(E + -E_ 




= (17) 

This assumption already entered the derivation of the system of wave equation, 
Eq. (|l0|), from the original Maxwell equations in the form of Eq. (|). However, 
this does not guarantee that all solutions of Eq. (|Io| ) or Eq. (|ll) satisfy Eq. 
([l7l). The latter is just the well-known statement that the electromagnetic field 
is purely transverse, ruling out longitudinal modes: the transverse electric field 
E± is related to the curl of the magnetic field by the Maxwell equation 

\7xB=-—=ikE ± , (18) 
c at 

and hence satisfies V ■ E± = 0; the longitudinal field Ei\, which can be written 
as the gradient of a potential <&, is responsible for violations of Eq. jl7| ). 

In view of the constraints imposed by the boundary conditions mm and by 
the zero divergence condition ( |l7j ) it is not possible (except for the case m = 0, 
as we shall see later) to set one (or two) of the vector components to zero without 
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setting the full electric field identically to zero, and thus it is not possible to 
reduce in a rigorous manner the vector problem into a scalar one. 

The problem can in principle be solved by converting Eq. ( |l7| ) from a condi- 
tion in the cavity volume to a boundary condition which can then be treated on 



the same footing as Eq. (16). One way of achieving this Q is by noting that if 

E = E ± +E {1 (19) 

fulfills Eq. (|l0|), then so does 

lE^VxVxSo^^VxVxii. (20) 

The latter is automatically divergence-free. In order for E to satisfy the bound- 
ary condition E t = 0, we require for Eq 

(e ) = andV -E Q = (21) 



t 

on the surface. Then one has indeed 



E t = r[VxVx E Q) 

= k 2 (^ +V(V- J Bo)) t = fc 2 (■$>) = 0- ( 22 ) 
on the boundary. The problem is therefore reduced to finding the auxiliary field 



E and then deducing the transverse field from Eq. (21). This leads to a system 
of three second-order differential equations for each vector component of Eq, 
given by Eq. (p"l|), all of which are coupled by boundary conditions that are, 
however, quite complex. 

The next step is then to write the field components as linear combinations 
of independent general solutions of Eq. ( |lTl) and determine the unknown coef- 
ficients in that expansion from the matching conditions at the boundary. The 
solution proceeds in an analogous but much less tedious way if we neglect the ad- 
ditional divergence condition. The important simplification is that we are then 
able to consider the E + and E_ components of the electric field independently, 
by setting all except for one component to zero. The boundary conditions ( |l6|) 
are then decoupled as well. More precisely, it will be shown that the wave equa- 
tions are then not only scalar but also separable, i.e., reducible to the solution 
of ordinary differential equations. 

We therefore would like to neglect the coupling that results from the condi- 
tion of zero divergence, provided that this can be justified in the context of the 
present study. There are various reasons why this approximation will provide 
us with useful results. Foremost, it will turn out below that the most important 
modes we find in this way in fact conspire to satisfy Eq. (ffTn a posteriori, cf. 



Section 3.5: the modes that provide the best confinement of the field in a tightly 



focused region around the focal point are the ones with m = 0. For these, the 
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different vector components decouple rigorously and the scalar program is exact. 
These m = modes are particularly significant because they provide the best 
confinement of the field in a tight focal volume. This is the paramount aim of 
the experimental dome structure. 

In addition to this exact result, the more transparent simplified problem 
allows us to evaluate the stability of the stationary states of the field in the 
parabolic cavity with respect to deviations from the confocality condition - a 
deformation that can readily occur in the course fabrication. This will be ad- 
dressed with the help of the ray picture in Section [?], and the ray trajectories 
themselves are independent of whether a vectorial or scalar field is considered. 
Since the exact nature of the deformation is unknown, it is necessary to make 
model assumptions and parametrize the deformation in some way. Although 
the range of possible behaviors explored within our model can be argued to 
be generic, we lose at that point the ability to predict accuratey all the indi- 
vidual modes of the specific sample. The error incurred by this fundamental 
uncertainty about the precise boundary shape is larger than the error made by 
adopting the simplified boundary conditions, and hence the latter are warranted 
on physical grounds. 

The consistency of these arguments is proven in Section ^ where we find 
that the only modes which can in fact be reliably predicted for a large range 
of possible deformations (because they are structurally stable against the emer- 
gence of chaos) are the ones with low m (or angular momentum in the classical 
picture), concentrated strongly near the z-axis. For these modes one can set 
approximately m » 0, E + = E- and E z = so that Eq. ( p~7| ) becomes valid. 

3.3 The wave equation in parabolic coordinates 

Having discussed the boundary conditions, we now provide the solutions to Eq. 
(^5|). In order to find a system of general solutions to the formally scalar differ- 
ential equations, Eqs. ( |i"l| c) or (|l5|), we express the scalar Laplacian appearing 
there in parabolic coordinates (£,77, </>), leading to the form 



dt; \ d£ J drj \ drj 



k 2 Q= — -Q (23) 
S?7 



where 



n = m + 1 for Q = E^ 



n = m-l for Q = E- (24) 
n = m for Q = E z 

Here we have used the fact that the derivative d 2 /d<j) 2 appearing in the Laplacian 
V 2 pulls down a factor —m 2 due to the ansatz Eq. ([12]). Although the righthand 
side is the analog of the centrifugal barrier in cylindrical problems, it thus 
depends not on angular momentum m directly but on a modified azimuthal 
mode number n. This occurs due to the additional d> derivatives introduced 
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when wc transformed the vector field components to cylindrical coordinates in 
Eq. (0). 

At this point we introduce the approximation of discarding the divergence 
condition so that we merely have to consider the boundary conditions ( |16j) w ith 
one and only one of the three field components nonzero. Then, Eq. (|23|) is 
separable in rj and £. We shall return to the details of the solution procedure 
in Section |J; for now it is sufficient to give the result: Denoting the separation 
constant by /3, the solution can be written in the form 

Q = F(k,0,S).F(k,-0,rj) (25) 

where F(k,(3,£) obeys 

ZF" + F'+(-^ + ¥-£ + ^F = (26) 

The functions F(k,f3,£) and F(k,—f3,rf) appearing here are solutions of this 
differential equation with the same k and n, but with sign-reversed (3, and hence 
their functional dependence on £ and r\ will be different unless (3 = 0. Without 
loss of generality we can assume n to be nonnegative, because it appears in the 
above equation only as n 2 . The solutions that do not diverge at £ = are of 
the form 

F(k, [3, £) = e^&M (^±1 -f;n+l; -ik^\ (27) 

where M(a, 6, z) is Kummer's confluent hypergeometric function. The function 
F as written here is in fact real, because of the Kummer transformation [fl3| 

M(^-a,b,-z) = e- z M(^ + a,b,z), (28) 

where we set a = i/3/k, b = n+ 1 and z = ik£. Appplying the theorem then 
yields F(k,p,0 =F(k,p,0*. 

The separation constant [3 and the wavenumber k at which to find the mode 
are still unknowns of the problem that have to be determined from the boundary 
conditions. The first constraint we can write down is 

F(fc,/3,2/) = (29) 

to enforce vanishing tangential field on the parabolic surface. In the two- 
dimensional plane spanned by the unknowns (3 and fc, this single equation defines 
a set of curves. The boundary condition on the focal plane requires that E± 
be symmetric under reflection, i.e., invariant under £ <-> r\. For E z , on the 
other hand, one needs odd parity. In order to construct such solutions with a 
well-defined parity, we have to form linear combinations 

E = F(k, & F(k, -13, n) ± F(k, -S3, F(k, (3, r,), (30) 
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where in addition 



F(fc,-/3,2/) = 0. 



(31) 



The set of curves parametrized by this constraint will intersect the curves defined 
by Eq. ( |2S| ) at certain isolated points in the (3 - k plane. By finding these 
intersection points, we determine the quantized values of (3 and k corresponding 
to solutions of Eq. ( p3| ) which satisfy the boundary conditions. It is not clear 
at this stage of the discussion how many intersections there are, or even how 
the curves defined by each equation separately will look. Before we analyze 
the different branches of these equations and identify their intersections based 
on asymptotic methods in Section |], it is useful to discuss in more detail the 
consistency of the fields thus obtained. 



3.4 Behavior at the focal point 

The main experimental purpose of the cavity is to concentrate the field near 
the focus as much as possible. Since one always has E z = there, it remains 
to discuss the behavior of E± in the focal region. Because of the "angular- 



momentum-barrier" on the righthand side of Eq. (23), the solutions F given in 
Eq. (|27| ) attain a factor (J 1 / 2 which suppresses the field near the origin £ = 
when n / 0. The Kummer function itself goes to M = 1 at £ = 0, so that 
the only way of getting a nonvanishing field at the origin is to set n — in 



Eq. (26). This means that the angular momentum quantum number must in 
fact satisfy m = 1 for E_ or m = — 1 for E + according to Eq. (§3). But this 
leads to a contradiction: if the field is nonzero at the origin, then because of 
the azimuthal factor exp(±i0) one faces a singularity at £ = r/ = in which the 
field is indeterminate. Therefore, there is no possibility to obtain a nonzero field 
precisely at the focus of the cavity. 

For m — ±1 there are still solutions of Eq. (E3), but they must involve 



solutions of Eq. (23) in suitable linear combinations such as to yield a vanishing 
field at £ = r) = 0. We have the freedom to linearly combine eigenstates of 
the wave equation at the same wavenumber k (yielding a stationary state with 
monochromatic time dependence) . First we use the real- valued solutions in Eq. 



( p7[ ) to form a superposition of the type Eq. (30) with a plus sign. Despite its 
symmetry it can also be made to vanish at £ = rj = 0, if one or both of the 
functions F(k, f3, £) and F(k, —(3, rf) are zero at the origin. 

Although we can hence find solutions for arbitrary m with a tangential elec- 
tric field that is symmetric under reflection at the focal plane, we can only 
attempt to concentrate the field near the focus, always with a node at the focal 
point, dictated for m ^ by the phase singularity at the origin. For m = 0, there 
is the residual angular momentum barrier due to n ^ 0, and thus even in this 
simple case - contrary to our expectation from quantum mechanical analogues 
- the "s-wave" solutions have vanishing field at the focus, as a consequence of 
the vector nature of the field. 
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3.5 Particular case: the fundamental s-wave 

The case m = can be discussed in more detail because it permits simple 
analytical expressions for the wave solution, if we specialize further to (3 = 0. 
In this case, the solutions in Eq. ( p7| ) simplify to 

F(fc,(U)cxI(!,^^"/ 2 , (32) 

dropping prefactors that are absorbed in the normalization. Here, I is the 
modified Bessel function. 

As was already noted below Eq. (|l5|), E± linearly dependent in the special 
case m = 0, so that we can in particular choose E + = E-. Then E p = and 
E + = E- = Ej,/V2- We thus obtain the TE field by setting 

E$ = Q = F(k,0,OHk,0,r)) (33) 

as in Eq. (p5|). This already satisfies the condition of symmetry with respect 
to the focal plane, without having to form a superposition of the type (|30|). 
Moreover, it satisfies the condition of vanishing divergence, as can be checked 
with Eq. |L7|). 

With n = 1 (for E + at m = 0), Eq. (|32"|) can be rewritten to obtain 




■ sin(fc£/2)sin(fc?7/2) 



(34) 



B 5 = -iE ■ ■ ^ ■ -±= sin(fc£/2) cos(fcr;/2) 
B = \ B, = +iE Q ■ ^ • J~L. . ^ cos(fce/2) sin(fc, ? /2) ( 35 ) 



and 



S = O 

The resonance condition is obtained from the boundary condition Eq. (j^) on 
the parabolic dome at £ = 2/ as 

k N = Nj, N -1,2,... (36) 

A relatively simple visualization of these modes can be obtained by express- 
ing the electric and magnetic fields in cylindrical and spherical coordinates which 
are more familiar. It can be verified using the relations between these coordi- 
nates to the parabolic variables that Eq. (34) then takes the form 

E p = 

E={ E z = (37) 
E 4> — E o ■ ^ ' (cos(fcz) - cos(fcr)) , 
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Figure 4: Cross-sectional view of the unfolded cavity with two closed, bowtic- 
shaped ray paths going through the common focus of the bounding parabolas. 
Families of such rays can be thought of as constituents of the s-wave in Eq. (|34]). 



and the corresponding magnetic field is 



B p = iE ■ ^(p • ^(zsin(fcr) - rsin(fcz)) 

B= { Bz =iE -^f s -^ (38) 
B^ = 

By splitting the various terms appearing here into two contributions, the elec- 
tromagnetic field can then be considered as the superposition of two fields: 

The first field is polarized along E$ and B p and can be expressed in cylin- 
drical coordinates as 

E E 
Ev~' — —5- cosik^z) and B^ = iJJlI—^— sin(/cjvz) (39) 
v knP ' k N p 

The second field is polarized in spherical coordinates (r, <f>, 9) along the direc- 
tions of the azimuthal and polar angles, <f> and 9, according to 

(2) _ E cos(k N r) (2) _ Eg sm(k N r) 

E * "sin0 k N r ^ B ° ~ Wfie sin9 k N r (4 ° } 

Here, we have used the substitution p — r sin# in the denominators. 

The first field, Eq. ([59|), corresponds to cylindrical standing waves with a 
phase variation along the z direction, while the second field, Eq. (f40|), cor- 
responds to spherical standing waves with a phase variation along the radial 
direction. This configuration is reminiscent of what is expected from a simple 
geometrical optics argument in which a ray bundle emerging from the focal point 
can propagate outwards as a spherical wave, upon reflection on the parabola it 
gets converted into a cylindrical wave, which in turn can counter-propagate back 
to the focal point after being reflected on the planar mirror and a second time 
on the parabola. In the unfolded double-paraboloid, the ray trajectories are of 
the type shown in Fig. ||. It should be noted that these two partial waves are 
not physical when taken individually, because in both cases the electric field 
diverges along the axis of the parabola. The divergence, however, cancels out 
when the superposition of the two partial waves is considered. We have so far 
only drawn this ray interpretation from a particular decomposition of the exact 
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solution; the question is how arbitrary that decomposition is and what we can 
learn from it. This will be pursued in the following section. The actual intensity 



distribution of these n = 1 states in the cavity will be plotted in Section 5.2 



where we can compare their spatial patterns with those obtained for larger n, 
in order to justify our claim that the s-wave modes provide the best focussing. 



4 Finding the modes within the short-wavelength 
approximation 

Having seen that even the long-wavenlcngth s-wave in our cavity can be inter- 
preted as standing waves arising from counterpropagating ray bundles and their 
accompanying wavefronts, we now turn to a more quantitative eikonal analysis. 
Such an analysis can provide accurate starting points for a numerical search of 
the exact wave solutions, which are determined by finding intersection points 
between the families of curves (|2^) and (|3l]) in the plane of (3 vs. k. Such 
semiclassical considerations based on the short-wavelength approximation are 
an important first step because there are, as we shall see, infinitely many in- 
tesections between the sets of curves determining the exact solutions, and one 
desires a means of finding them in a systematic way, labeling them by "quantum 
numbers" , giving the number of nodes in the field along the coordinate lines for 
£ and r\. Beyond this very practical use of the short-wavelength limit, we also 
want to establish a physical understanding of the resonator modes that allows 
us to predict how they depend on changes in the cavity shape. This aspect of 
the ray picture will be expounded in the last section. 



4.1 WKB approximation and effective potential 



The equation to be solved is Eq. (|26|), an ordinary second-order differential 
equation, where the angular momentum m enters as a parameter trough the 
constant n. We are looking for solutions F(£) which satisfy the boundary con- 
dition F(2 f) = and are not singular at £ = 0. The standard short-wavelength 
approach to be employed here is the WKB approximation jT4|. After division 



by £, Eq. (26) takes the form 

d 2 f ldf 1 A 2 4/3 n 2 \ 

^ + ^ + + T'e) m = - (41) 

For the subsequent analysis it is convenient to introduce a dimensionless coor- 
dinate 

u = y/kl. (42) 

Dividing Eq. ( |4l] ) by k 2 and defining a rescaled separation constant 

(43) 
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Figure 5: Solid curve: the effective potential V(u) for n = 3, showing the 
classical turning points uq (vq) where Z (-Z) intersects the effective potential 
V, cf. the dashed lines. The ranges of classically allowed motion for the two 
degrees of freedom u and v with energies Z and —Z are indicated by the shaded 
bars (dark for u, light for v). The outer turning points at u, v = \J2kj act as a 
hard wall whose position depends on k. 

the following equation is obtained: 

d 2 f 1 df 1 fn 2 a \ . rr ... 

This has a form similar to the one-dimensional Schrodinger equation of quan- 
tum mechanics, except for the first u-derivative which makes the kinetic energy 
operator non-selfadjoint. This term appears in the radial equation of cylin- 
drically symmetric problems but does not affect the applicability of the WKB 
approximation p5| . 

The WKB solution requires us to find the classical turning points in the 
potential appearing in this equation, with Z playing the role of the total energy. 
This effective potential, 

V(u) = l(^-u*y (45) 

is a superposition of an inverted parabola and the centrifugal potential deter- 
mined by n, giving rise to the solid line in Fig. |[ Using this together with the 
ansatz 

/(„)„ 1 P (u) = v /Z-V(u), (46) 

p(u) 

the approximate solutions are found by imposing the boundary conditions at 
the turning points. 

There is only one possible turning point corresponding to the closest ap- 
proach to the origin u = 0, which is given by 



V(u )=Z^u = ^n 2 + ^-^. (47) 

If Z > and n = 0, then no inner turning point exists. This inner turning point, 
in classical mechanics, is the point where the momentum in the x-direction 
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smoothly goes through zero as it changes sign, and hence the probability per 
unit time ol finding the particle becomes infinite. In the ray dynamics, this 
phenomenon gives rise to a caustic. This will be discussed further in Section ^. 

The outer turning point of this classical picture is determined by the Dirichlet 
boundary condition at the parabolic mirror, which in the new coordinate is 
located at 

& = 2/=> Ul = y/2kf. (48) 



It is the boundary condition f{u\) = f(y/2 kf) = in which the short-wavelength 
condition is contained: we assume that at the outer boundary the wavefunction 
has the WKB form, Eq. ([Iq), which requires that the dimensionless size param- 
eter satisfies 

x=2fc/»l, (49) 

i.e. this additional boundary is far away from the classical turning point uq of 
the effective potential. All steps discussed above for /(£) apply analogously to 
the variable rj appearing in the product ansatz Q, Eq. (p5|), if we reverse the 
sign of Z everywhere and replace u by the variable 

v = y/kq. (50) 

Then the inner turning point vq for this second degree of freedom is obtained 

as 



V{v ) = -Z^v = ^v? + — + -. (51) 
The values of uq and vq determine the distance of closest approach to the z axis 




4.2 Quantization conditions 

Under this condition, the semiclassical determination of the eigenfrequencies 
proceeds by applying the Bohr-Sommerfeld quantization to the action integral 
for one period of the motion in the effective potential. One round trip consists 
of the path from uq to iti and back to uq. The quantized action is therefore 

J(Z, x;n,v) = 2 J y/z - V(u) du = 2tt (v + |j . (52) 

The integer v = 0, 1 . . . is the number of nodes of the wavefunction in the 
potential, and the constant 3/4 takes into account the phase shifts of 7r and 
7r/2, at the outer and inner turning points, respectively. In other words, the 
above quantization condition is an approximate way of writing the phase-shift 
requirements that hold at boundaries and caustics, using the approximation that 
the wave propagation itself is described by a wavefront whose phase advance in 
x is given by the function J. 
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The result of the integration is found to be 
J(Z, x; n, v) — \J x 2 + Z x — n 2 



Z 1 Vx 2 + Zx-n 2 + X + Z/2 
2 11 y/n 2 + Z 2 /A 

Zx-2n 2 7T 



n arcsm ■ 



^Z 2 + An 2 2 



x 

- 2tt (v+iy (53) 

This is an equation for the two unknowns Z and x, i.e. for the rescaled separation 
constant and size parameter. It is therefore analogous to Eq. (|2^). The integer 
v uniquely labels all the allowed solutions /„ of Eq. (Q). This is an important 
difference to Eq. (p9|): there, the function F in fact has infinitely many branches 
that satisfy the equation, which are however not labeled explicitly. The great 
advantage of Eq. ( |53| ) is that these branches are explicitly enumerated by so 
that fixing this index selects exactly one curve in the Z-x plane instead of an 
infinte family. 

As in Eq. ( |25| ) , the field consists of products of the form 

Q = f u (x,Z,u)f fi (x,-Z,v) (54) 

with v defined as in Eq. (|5^). The two function /„ and / M have their analog in 
the exact solutions F of Eq. (|26|), corresponding to the branches of F labeled 
by v and /i, respectively. 

We then form combinations the form of Eq. (|3(]) to enforce the required 
symmetry with respect to the focal plane. The semiclassical WKB quantization 
for the function /^(x, —Z, v) provides a second equation of the form (|5^), 



J(-Z,x;n,fi) = 2tt I fi+ - J . (55) 

These two quantization conditions play the same role as Eqs. (||) and (|l]): the 
intersections of the curves parametrized by them determine the quantized values 
of Z and x. However, the WKB method affords a great simplification: by fixing 
the branches v and /x, the intersection of the two resulting curves is uniquely 
determined. To illustrate this situation, we show in Fig. || how the lines defined 
by the above two equations traverse the Z-x plane. Only a small portion of this 
plane is shown, emphasizing the behavior of the semiclassical results at small 
x where their accuracy should be at a minimum. Comparison with the exact 
families of curves shows, however, that the WKB results are excellent even in this 
long- wavelength limit. Note that by symmetry, intersections occurring at Z = 
are always between curves with the same branch index v = ji. All curves shift 
to larger x with increasing n because of the larger centrifugal barrier, pushing 
the classically allowed regions of the effective potential in Fig. [sj outward. 
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Figure 6: For the first six values of n, the graphical solution of the simultaneous 
equations ( |53|) and ([55]) can be read off from the intersections of the red curves. 
The dashed black curves show the analogous graphical solution of Eqs. ( ^9|) and 
(f3ll). The exact and semiclassical curves are almost indistinguishable (except 
for n = 0), attesting to the striking accuracy of the former even at the smallest 
possible size parameters x. All plots can be continued to Z > by reflecting at 
the axis Z = 0. The WKB curves with positive slope belong to Eq. (1551), the 
falling lines are created by Eq. (^3|). They are labeled starting from the leftmost 
by fi, v = 0, 1, 2. ... 
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Figure 7: Four states with n = 1. Grayscale indicates the magnitude of the 
electric field for the modes (E±) of the confocal double paraboloid, highest fields 
shown in black. The vertical axis is z, the horizontal axis the axial distance p. 
Increasing x means shorter wavelength and hence more nodal lines (white). The 
size parameter is quantized according to Eq. (^6|) with N = 1, 2 in the top row, 
and N = 7, 14 at the bottom. 



5 Exact solution for the modes and their field 
distribution 

Once the allowed combinations of Z and x - or equivalently (3 and k - at which 
the boundary conditions are satisfied have been found, the problem of finding 
the modes is solved. For example, we can now plot the intensity distribution of 
each mode by using the quantized values of (3 and k in Eq. ( |2~^ ) and forming 
the proper linear combinations of the form Eq. (|3^) . 



5.1 Mode profiles 

This will now be carried out for the lowest-lying modes as obtained from the 
intersections in Fig. [| Any given value of n can have a different meaning for 
the intensity distribution in the azimuthal direction, depending on which case 
in Eq. ( |24| ) we choose to consider: n = m± 1 for the modes. Since the azimuthal 
field variation is trivial, tx exp(zm0), we wish to restrict our attention to the 
mode profile in the plane spanned by p and z in cylindrical coordinates. The 
variable governing this property is n, not m. Therefore, n is used here to classify 
the mode profiles. 

As has been done in the previous sections, we shall take the focal plane to 
be the symmetry plane of a double paraboloid, and plot the wave fields in this 
unfolded cavity. This is done in view of the subsequent discussion, where we 
shall establish the connection between these modes and the ray dynamics. Some 
wave plots are shown in Figs. and ||. Note that the case n = does not appear 
among the solutions listed here because it corresponds to wave fields that do 
not vanish on the z-axis and hence are irreconcilable with the finite angular 



momentum m = ±1, as discussed in section 3.4 ) 



If we look at only the left column of Figs. ^ an d tl h is apparent that all 
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Figure 8: The modes shown here do not follow the simple law of Eq. ( pq ) but 
were obtained numerically. With increasing centrifugal barrier, labeled by n, 
the forbidden region around the z-axis grows outward. 



states with Z = look similar, as do all states with Z ^ 0. A similar observation 
can be made in the right columns of the figures. Comparison to the intersecting 
lines in the graphical solution, Fig. [| shows that states with the same nodal 
pattern indeed result from the crossing of the same pair of lines - labeled by 
the same /i and v 1 only for different n which pushes the intersecting lines to 
higher x. However, the interpretation of fi, v as the number of nodes in the 
parabolic coordinate directions cannot b e ca rried through in all of the plots. 
We will return to this problem in Section 5.3. 

In order to achieve the best possible concentration of fields near the focus, 
the most promising candidates are the modes with n — 1. Among these, the 
patterns shown in Fig. fj] indicate that the states at Z = in turn show the 
highest intensity near the focal plane. These are precisely the fundamental s- 
waves we discussed in Section 3J5, with the wavenumbers quantized according 
to Eq. (pq), which for the size parameter reads 



x N = 2k,J = 2ttN. 



(56) 



This is an exact result which can be compared to the WKB quantization condi- 
tion in Eq. j5^ ) with Z = 0, n = 1. The latter actually has a more complicated 
form, 

\fx 2 - 1 + arcsin - = 2tt (v + 1), (57) 
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but to second order in the small quantity 1/x this is identical to Eq. ( p6| ) with 
N = v + 1. This confirms the observation made in Fig. || that the numeri- 
cal agreement between exact and semiclassical solutions is good even for small 
quantum numbers. 



5.2 Focussing and the effective mode volume 

In order to evaluate the field enhancement that is achieved in the fundamental 
TE modes discussed in Section 3.5 and shown in Fig. |?], it is necessary to examine 



the distribution of the electromagnetic energy in that mode. The energy in a 
parabolic cavity of focal length / is 



U 



H=2f rV =2f r 4>=2-K 



4 j£ =0 J v =o Jcfr=o 
|(|^| 2 + \E V \ 2 + \E^\ 2 ) + i-(|5 € | 2 + \B n \ 2 + 

^p-d^drjdt/) (58) 



which gives, for the fundamental (s-wave) TE modes, Eq. (p4[), 

where the value of the integral can be evaluated numerically. 

For the experimentally realized cavity described in Section 1, kf = 14 it 
so that the value of the integral is 2.527. The intensity distribution for this 
mode is shown at the bottom right of Fig. (note x — 2kf). To examine the 
energy distribution in the cavity, we can evaluate the energy that is contained 
at each lobe of the standing wave of parabolic wavefront that corresponds to the 
mode. We note then that the first lobe, corresponding to a parabolic wavefront 
of focal length /i = A/2, contains 48% of the total energy; to see this, replace 
the integration limit in Eq. ( |59| ) by n. This lobe occupies a physical volume of 
Vq = 7rA 3 /4 whereas the volume of the overall cavity, V = 2ir f 3 , is 2744 times 
larger. 

This underscores the very large confinement of the field that occurs in the 
vicinity of the focal point and points to the possibility of observing a very large 
enhancement of spontaneous emission into this mode. The fraction of the total 
energy contained in the first lobe of course reaches 100% if the smallest possible 
cavity with k f — n is considered. However, the size achieved in our present 
sample already approaches the optimal conditions if one takes into account that 
enhancement of spontaneous emission requires not only a small effective mode 
volume but most of all a small local density of states | ft6| . The average density 
of modes in an arbitrarily-shaped electromagnetic resonator of volume V is a 
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fundamental quantity in the theory of blackbody radiation and was derived by 



indicating that the number of modes in the interval dk is proportional to the 
number of additional volume quanta Vq that fit into the given volume V when 
k increases to k + dk. The local density of states in the focal volume Vq can 
therefore be interpreted to be the same as the total density of states in a small 
cavity of volume V = Vq. This, in turn, is roughly the effective mode volume 
for the fundamental s-wave in our structure. /From this we conclude that the 
spontaneous emission enhancement should be close to the maximum possible 
value even though our cavity is not of the minimum possible size. This is one 
of the central advantages we were looking for in the parabolic cavity design. In 
this discussion we have assumed for simplicity that the Q-factor of the modes 
under consideration is fixed, independent of size and quantum numbers. This 
severe simplification will be removed in Section |^. 

In the higher order modes with m > 0, the centrifugal barrier prevents the 
field from approaching the focal point. This implies that these modes will 
have a larger effective volume and, correspondingly, a smaller enhancement 
of the spontaneous emission rate. An added difficulty concerning the higher 
order modes arises from the limited experimental control over the exact cavity 
shape. As discussed in Section [?], small deformations of the cavity (modeled as 
deviations from confocality) , result in chaos, leading to a loss of constraints on 
the possible regions of phase space which can be explored. This further increases 
the effective volume of these modes. The enhanced spontaneous emission into 
the fundamental s-wave implies that this mode will also exhibit a large gain 
and, correspondingly, a low lasing threshold. The preliminary conclusion of this 
section is therefore that a mode with low angular momentum and small Z (or 
(3) will be the dominant mode in a laser of parabolic geometry. 

5.3 Caustic structure in the wave solutions 

In order to arrive at the solutions shown in Fig. |^, we started from the semi- 
classical (short-wavelength) approximation and then refined the quantized Z 
and k further by applying the exact modal conditions. However, the question 
arises how the quantum numbers fi and v which label the semiclassical solutions 
can be visualized in Fig. |8[ The answer is that the symmetrization procedure 
obscures this identification. What happens can be understood if we ignore the 
parity requirement and plot the wave fields in the simple product form of Eq. 




(60) 




(61) 
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Figure 9: Mode intensities as in Fig. |8j, but without performing the symmetriza- 
tion prescribed by the focal plane boundary condition. The WKB quantum 
numbers \i (and v) can be read off by counting the number of wavefunction 
nodes parallel (and perpendicular) to the reference line 7. Modes in the first 
three rows correspond to the symmetrized versions of Fig. ||. In order to illus- 
trate the approach to the short- wavelength limit, additional modes are shown for 
which the formation of caustics is apparent in the high-intensity ridges (black) 
bordering the classically forbidden regions (white) . 

The symmetrization performed according to Eq. ( |30|) with A = B intro- 
duces no change whatsoever if the separation constant is /3 = 0. Therefore, the 
intensity profiles of all modes with Z — in Fig. || are the same before and 
after symmetrization. However, the wave patterns acquire a qualitatively differ- 
ent and simpler form if we desymmetrize the remaining states. This is shown in 
Fig. ^| The nodal patterns now appear in a regular fashion along the coordinate 
lines for u and v (or £ and 77), and their number along these lines is uniquely 
determined by fi and v. 

By symmetrizing a state such as the one shown for n = 3, k = 18.59 and 
Z = —2.54 in Fig. |9[ the field shown in the desymmetrized plot is added to 
its reflection at the focal plane, thus allowing some nodal lines to be "filled 
in" , as seen in the corresponding state at the bottom right of Fig. ||. The 
desymmetrized waves in Fig. ^ exhibit nodal lines precisely along lines of 77 = 
const or £ = const. In addition to the simple nodal structure, we also observe a 
clear segregation between regions of negligible intensity and regions of oscillatory 
field, with dividing lines between them that become more and more pronounced 
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as the size parameter x — 2 k f increases. These are the caustics, which in fact 
accumulate an increasing amount of intensity as the short-wavelength limit is 
approached. The caustics follow parabolic coordinate lines as well, as is apparent 
from the last row of Fig. |. The field at n = 3, x = 34.16, Z = -8.72 is bounded 
from below by a broad inverted parabola, and excluded from the z-axis by a 
steep upright parabola. The intersection of both parabolas forms the caustic. In 
the mode at n — 4, x — 73.73, Z — 0, both the upright and inverted bounding 
parabolas are symmetric as we expect for Z = 0. 



6 Caustic structure in the ray picture 

In this section we will elaborate on the relation between mode structure and ray 
dynamics, as a basis on which we can predict the effect of shape perturbations 
on the mode structure. The caustic patterns revealed in the last section by the 
decomposition into the product states as in Eq. ( |25| ) is a direct consequence of 
the classical turning points in the effective potential V, Eq. (ff5|), for the motion 
along the £ and r\ directions. The distinction between classically allowed and 
forbidden regions gives rise to the regions of oscillatory and vanishing fields in 
Fig. H The effective potential has, so far, been discussed only as an auxiliary 
concept that proved convenient in the WKB treatment; its relation to the be- 
havior of the rays of geometric optics is, however, well-known. For the sake of a 
self-contained presentation, we convey here the idea behind the general eikonal 
theory by showing how to derive ray equations from the one-dimensional sepa- 



rated wave equations, Eq. (26). The argument is non-standard in the sense that 
Eq. ( p6f ) is based on the full vectorial wave equation (i.e. with polarization), 
and we therefore shall find that for a given angular momentum m, slightly dif- 
ferent ray trajectories have to be considered depending on polarization. This is 
because the quantity entering Eq. (Bq) is n, not m. 



6.1 Ray equations from the WKB approximation 

Inserting the WKB ansatz, Eq. (ff6|), into the wave equation for the separated 
variables, Eq. (||^), one finds to leading semiclassical order that p must satisfy 
the equation 

p 2 u + V(u) = Z, similarly p 2 v + V(v) = -Z. (62) 

We can interpret this as the Hamiltonians of two decoupled linear systems, and 
add them to obtain the Hamiltonian for the combined system, 

H =pl+p 2 v + V(u) + V(v) (63) 

The trajectories we are looking for then satisfy the equation H(p u ,p v ,u,v) = 
Z — Z = 0, or written out: 

1 / n 2 ?i 2 \ 
Pl+Pl + l {^ + ^- u2 -v 2 ) =0. (64) 
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If we divide this by (u 2 + v 2 ), the result is 



p2u+T>2v 1 ' " 1 ) = 0. (65) 



u 2 + v 2 4 

this can also be interpreted as arising from a new Hamiltonian 

u z + v z 4 u v 

by requiring 

H(p u ,Pv,u,v) = -. (67) 

The Hamiltonian in this form is analogous to the wave equation in parabolic 
coordinates, Eq. (p3|), where the Laplacian is divided by the same scale factor 
(u 2 + v 2 ) that accompanies the conjugate momenta here. One can now use 
Hamilton's equation of motion to replace momenta by "velocities" , the definition 
being 

dH dH 



OPu OPv 



(68) 



This leads to the substitution 



Pu = ^ (u 2 + v 2 ) it, p v = i (u 2 + v 2 ) u, (69) 
which brings Eq. (|67| ) into the form 

(u 2 + v 2 )(u 2 +v 2 ) + ^ = l. (70) 
Reverting to cylinder coordinates, the above equation becomes 

k 2 ( P 2 + * 2 ) + £- p = l. (71) 

Here we used the definitions of the coordinates in Eqs. (||) and (f42|). To examine 
what this equation has to do with the ray dynamics, we take the ray-picture 
point of view now. 



6.2 Geometric optics in cylindrical coordinates 

If we consider the three-dimensional motion of rays in a double paraboloid of the 
shape in Fig. || (b), their propagation between reflections at the parabolic walls 
will of course follow straight lines, and hence there is no place for any coordinate- 
dependent potential V. However, in order to compare the ray dynamics to 
wavefunction plots in the z-p plane as shown in Fig. ^, we must project the 
ray motion onto this plane as well. In the wave analysis, this projection was 
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achieved by using the rotational symmetry of the cavity around the z axis to 
eliminate the azimuthal coordinate <fi from the problem in favor of the angular 
momentum quantum number m. 

In ray optics, we can do the same: rays can be classified by an angular 
momentum L z because of the axial symmetry. To see this, we first define L z . 
A ray trajectory is a curve consisting of straight line segments between each 
reflection. If we parametrize this curve as r(Z), where I is the path length along 
the ray from some arbitrary starting point, then \r(l)\ = 1. Here and in the 
following, the dot represents the differentiation with respect to arc length, d/dl. 
In cylinder coordinates p, <f>, z, we can decompose this as 

r = pe p + ze z + r 4>e^. (72) 

Between any two reflections, this is a constant unit vector in the direction of 
the ray. With this, the equation for a straight line segment can be written in 
general as 

r x r = L, (73) 

where L is a constant analogous to the angular momentum of classical mechan- 
ics. 

Because of the rotational symmetry around the z-axis, the azimuthal unit 
vector at the point of reflection is always tangent to the surface. Therefore, a 
reflection does not change the component of r along e^, so that p4> is continuous. 
Since the ray curve is itself continuous everywhere, so is p(l). Hence the quantity 

L z = p 2 i> (74) 

is also continuous at each reflection. But this is just the z-component of L in 
Eq. (|73|), as can be verified by performing the cross product there. Thus, L z is a 
constant between reflections, which together with its overall continuity implies 
that it is a conserved quantity for the whole ray trajectory. 

Using Eq. (f74|), the fact that r is a unit vector, Eq. (|72|), can be recast as 

P 2 + z 2 + % = 1. (75) 
P 2 

^,From the ray approach we have thus obtained an equation almost identical 
to Eq. (^Tj). We only have to re-define the path length variable / to make it 
dimensionless, by introducing 

s = kl, (76) 



to obtain for the derivatives 



and interpret furthermore 



L z = l (78) 
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Then Eqs. ( |7l| ) and (|7q) become identical, if we interpret the dot in Eq. ([71]) 
to mean d/ds. The scale factor of the "time" variable parametrizing our trajec- 
tories is irrelevant for the shape of the paths, so that we can conclude that the 
ray picture introduced here is equivalent to the motion described by the WKB 
effective potential, with the important identification of Eq. ft7i\). 

Besides Eq. ((75|), the only other equation that is needed to completely de- 
termine any ray trajectory from its initial conditions is the law of specular 
reflection, which can be formulated with the help of the outward normal unit 
vector u at the reflection point as 

Yre fleeted = T - 2 U (u • f). (79) 

This corresponds to a reversal of the normal component of r. Here we can 
see explicitly that reflections do not affect the component of r in the direction 
of e^, since the normal u has no e^-component as a consequence of the axial 
symmetry. 

This latter fact also means that we can simply drop the e^-component from 
Eq. (f79]) altogether. Therefore, we now define the two-component vectors in 
the z - p plane by dropping the e^-components from the corresponding three- 
component vectors. Thus, r becomes 

v = />e, + ie,= fP J , (80) 



and similarly 



u = Up e p + u z e z . (81) 



In this two-dimensional space, the specular-reflection condition retains the form 
ofEq. ©, 

^reflected = V - 2 U (u • v). (82) 

This is the reason why we can call the motion in the z - p plane a billiard 
problem. 

6.3 Curved ray paths in the centrifugal billiard 

We know that the trajectories between reflections are straight lines, so that the 
components of r in the cartesian coordinate frame are constant for each segment. 
In our new z - p frame of reference, the z-axis is the same as the cartesian one, 
so that we still have v z = i = const between reflections in Eq. (|so|). However, 
the same does not hold for the p-component of v. Instead, we obtain from Eq. 
75b 

L 2 

p 2 + -§■ = I - z 2 = const. (83) 
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If we multiply this by 4p 2 , it can be written as a differential equation for p 2 : 

{j/} ! = W = 4(l-iV-4£*. (84) 

The solution is that p 2 (l) describes a shifted parabola, 

p 2 {l) = pi + 2 (I - h) ^f(l-z 2 )p 2 -L 2 + (l - k) 2 (1 - i 2 ), (85) 

where p 2 is the integration constant and specifies the value of p 2 (h) at the 
starting point h of the ray. Since furthermore z is a linear function of I (i = 
const), we can for definiteness fix the initial point is to lie on the focal plane 
and substitute 

I = l(z) = z/z (86) 

to find that Eq. ( |85| ) describes a curved path p(z) in the z -p plane. The curved 
nature of this trajectory is a direct consequence of the centrifugal potential 
L 2 /p 2 in Eq. (|75j), and we would recover straight lines, i.e. linear variation of 
p(l) for L z = 0. This is why we refer to this problem as a centrifugal billiard |fl7| . 
For a visual example of how curved traces arise from stright-line trajectories, 
the reader is referred to Fig. |l3| (e) which will be discussed in Section |t]. 

An example of the ray motion in the special case L z = is already displayed 
in Fig. |], showing no curved trajectories because there the z - p plane is indis- 
tinguishable from the cartesian z - x plane. For L z ^ 0, curved ray trajectories 
in the z - p plane are shown in Fig. [l^ for four different initial conditions under 
which the ray is launched. Note that the parameter L z as given in the plot 
has dimensions of length, cf. Eq. (f75|). This reduced two-dimensional problem 
can be analyzed completely without reference to the original three-dimensional 
ray tracing, with L z as a parameter that encapsulates the third degree of free- 
dom <j> which has been eliminated. We have, broadly speaking, converted to a 
co-rotating frame of reference (with rotation speed always matching the vary- 
ing angular velocity of the ray), and thus obtained a planar problem in which 
we now look for the classical orbits. The simplification is considerable because 
the three-dimensional ray motion in the cavity is rather difficult to visualize, 
compared to the motion in the z - p plane. 

The two periodic orbits in Fig. [n] (a) exemplify this situation: after the 
trajectory completes one round-trip in the z - p plane, it returns to its initial 
position with the initial orientation - but in the original three-dimensional cavity 
there has also been a motion in <j> which does not necessarily amount to a full 
rotation around the z axis. Hence, this periodic orbit of the centrifugal billiard 
is not in general a true periodic orbit of the parabolic dome, cf. Fig. ^ (e). 
However, we can reverse this statement and conclude that any periodic orbit of 
the three-dimensional problem must also be periodic in the z - p plane. This 
cautionary remark concerning the interpretation of Fig. ^ is relevant if we 
attempt to interpret the actual modes of the original cavity in terms of a naive 
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Figure 10: All starting conditions for the rays in (a), (b) and (c) amount to the 
same angular momentum, L z = 0.1 /, but prescribe different angles of incidence 
with respect to the boundary, and different positions of impact, (a) shows 
two different orbits in the same plot, both are periodic and symmetric in this 
projection onto the z - p plane, differing only in their starting points, (b) 
shows a single path, which is quasi-periodic because it does not close on itself 
even in this projection. Instead, it densely fills a region of space delimited by 
a caustic whose shape is given by the parabolas (blue). The caustic becomes 
more asymmetric in (c) and (d), each of which shows a quasiperiodic orbit. 
The angular momentum in (d) is L z — 0.6/, leading to a larger forbidden zone 
around the z - axis. 



physical optics approach: one might think that a quantized mode has to be 
associated with ray paths that form a closed loop and in that way "reproduce" 
themselves. However, a comparison between Figs. [| and [l^ reveals that periodic 
orbits seem to play no special role for the mode structure. 

What shapes the modes is not any single periodic ray orbit, but the caustics 
as they appear in Fig. [l(] (b) - (d). The spatial distribution of the ray trajectories 
exhibits a clear correspondence with the modal intensities shown in Fig. ^|, 
particularly in the shape of the caustics. This is most convincing for the two 
examples in the bottom row of Fig. ^ where the wavelength is shortest: The 
state at n = 3, x — 34.16, Z = —8.72 should be compared to Fig. [lO] (c), and the 
reflection-symmetric mode with n = 4, x — 73.73, Z = finds its counterpart 
in Fig. 0(d). 

Caustics are immediately generated when we follow a single quasiperiodic 
orbit, but not so for a periodic one. However, periodic orbits occur in infintc 
families which, when plotted together, again fill a region of space bounded by a 
caustic curve. The two members of the family shown in Fig. [l^ (a) are obtained 
by launching a ray from the focal plane, perpendicular to it, differing only in 



the radial distance p of the launch. All other siblings of the examples in Fig. 10 
(a) combined, would create a picture almost identical to the one generated by 
the single quasiperiodic orbit in Fig. [lO] (b) - the latter is in fact the result of 
only a slight deviation from the initial conditions chosen in Fig. ^ (a), with the 
result that the orbit almost, but not quite, closes on itself after one round trip, 
and continues to fall short of closing itself after each subsequent round trip as 
well. The conclusion is that from the point of view of the caustic structure in 
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our system, there is no qualitative difference between periodic and quasiperiodic 
orbits. 

The fact that all orbits can be characterized by a particular caustic which 
they touch, and that moreover all periodic orbits come in infinite families, is a 
general property of integrable Hamiltonian systems, to which the special cen- 
trifugal billiard defined here belongs. That the paraboloid billiard is integrable, 
can already be concluded from the existence of a separation ansatz for the wave 
equation, which we discussed in Section ||. However, we have not yet completed 
our program of connecting the ray and wave approaches, and in particular we 
have not addressed the question of how to determine quantitatively the type of 
ray trajectories that correspond to a given mode. So far, the correspondence 
was established by visual inspection alone. The quantitative connection is ob- 
tained by comparing the ray patterns of Fig. [ijj with the effective potential V 
of Eq. (f45|). We shall see that for an integrable system, we can in fact uniquely 
connect a particular caustic with a given mode. 

As a final remark concerning the periodic orbits in this integrable system, 
it is worth comparing the patterns of Fig. HQ (a) and especially Fig. (0) with 



the "bowtie laser" of Ref. |18 . There, a semiconductor cavity was designed in 
such a way as to obtain lasing from a bowtie-shaped mode with highly desirable 
properties, foremost among them its focussing action in the center of the cavity. 
The focussing patterns of Fig. ([!]) are very similar, but the main difference is 
that in our case these orbits occur in families whose members can cross the z = 
plane with all possible axial displacements p. In the semiconductor cavity, most 
rays move on chaotic trajectories, and only a small range of initial conditions 
for the rays lead to a stable bowtie pattern, leading to modes which are strongly 
concentrated near a unique bowtie path, and hence even less spread out in space 
than the examples shown in Fig. [7|. This leads us to anticipate that the beneficial 
properties of the n = 1 modes found for our integrable system can in fact be 
enhanced if we allow for the possibility of chaos in the ray dynamics. 

6.4 Connection with the effective potential in parabolic 
coordinates 

The classical turning points for the two degrees of freedom u and v in the 
potential of Eq. (^|) determine the parabolas which describe the caustics in Fig. 
O . We notice that the caustics (and also the quasiperiodic rays that generate 
the caustics we show) have a well-defined distance of closest approach po with 
respect to the z axis, given by the corner at which the two bounding parabolas 
meet. Describing this in parabolic coordinates, we find that po is approached if 
both £ and r\ simultaneously reach their inner turning points. Expressing this 
condition in terms of Eqs. ( [f7| ) and (jsil), we obtain the simple semiclassical 
relation 

(87) 
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Here, we have used the coordinate transformation p = y/£ 7], cf. Eq. (^), and the 
definition of the rescaled variables, Eq. ji^). 

The distance of closest approach for individual periodic orbits is not given 
by this expression, but the minimal p over the whole family of such orbits does 
follow this law. The caustics in Figs. |l0| (b-d) exhibit cusp singularities at 
po because in that extreme point the z motion has zero velocity: it is clear 
from Eq. j7^ ) that the smallest p will be achieved when p = z = 0. But 
from the same equation we immediately obtain that the angular momentum 
then equals the axial distance, and with Eq. ( |87| ) this reproduces Eq. (|78|). 
We have therefore established that the ray's "angular momentum" is directly 
proportional to the modified angular momentum quantum number n of the 
mode under consideration. In the semiclassical limit of large k, the difference 
between n = m ± 1 and m becomes negligible in this expression, so that we 
recover the intuitively expected proportionality 

777 

L z = p ~j. (88) 

This approximation means that we can neglect the effect of polarization on the 
ray-wave correspondence in the semiclassical limit - however, we shall make use 
of this only later, in the ray analysis of Section ^. Since we have been interested 



in states at rather small k and in particular n > 1, we have plotted in Fig. 10 
only trajectories with L z ^ 



A second semiclassical relation follows from Eqs. ( |47| ) and (51) if we ask for 
the value zq of z corresponding to the point Uo, vq at which the caustics have 
their singularities. The whole caustic is uniquely determined by its singular 
point at radial distance p and height z , cf. Fig. [l0| According to Eq. (||), we 
get 

* = f . (89) 

This identifies the meaning of the separation constant Z, also quantifying the 
earlier observation that for E = both the wave and ray patterns are symmetric 
with respect to the focal plane: in that case, the cusp occurs on this mirror plane, 
as in Fig. |l^ (b) 



With Eqs. (|78|) and (89), we have completed the bridge from the exact wave 
equation via semiclassical WKB quantization to the ray caustics. By specifying 
the quantized n, Z and A: of a given mode, we uniquely determine a caustic and 
with it a particular family of ray paths. Now we can use additional properties 
of the ray picture to better understand the cavity modes. This is especially 
promising in this system because we have seen that the semiclassical approx- 
imation is extremely accurate here. The reason for this somewhat surprising 
accuracy lies itself in the properties of the ray dynamics, but in order to make 
this clearer we need to introduce the concept of a phase space in which the ray 
dynamics can be described. 
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6.5 Families of rays and Poincare sections 

A phase-space description is often used in classical mechanics because it carries 
more information about the possible trajectories than mere real-space diagrams. 
This approach has recently been applied to the analysis of ray dynamics in 



optical cavities as well 26, [27|], with the goal of providing insights that are 
not revealed by ray tracing in real space. In particular for the treatment of 
non-integrable resonator geometries, it has proved valuable to represent the 
phase space of the rays in terms of Poincare Surfaces of Section (SOS). For our 
purposes, the following SOS will be chosen: 

It is easy to convince ourselves by recalling Fig. ^ that any ray trajectory in 
the cavity has to encounter the focal plane infinitely many times as it propagates. 
However, the radial distance of these crossings, as well as the value of p may 
vary from one crossing of this plane to the next. Now we can consider 

p and p p = p (90) 

as a pair of canonically conjugate position and momentum variables, and at- 
tempt to image the subset of phase space spanned by them. In order to do 
that, we launch a ray trajectory and follow it for many crossings of z — 0, each 
time recording the instantaneous values of p, p p as a point in a two-dimensional 



graph. The result is shown in Fig. 11. A typical trajectory is - as mentioned 
above - quasiperiodic, and in the SOS generates a dense set of points that all 
lie on a smooth curve. Several trajectories have been followed in this way and 
are represented in Fig. |ll| by the different individual curves. Each curve ex- 
hibits some minimal axial distance p m i n > po] this is a true inequality because 
quasiperiodic orbits do not reach their point of closest approach to the z axis 
precisely on the focal plane. Since the SOS records the instantaneous p upon 
crossing the focal plane, the resulting curves have their turning points at larger 
P- 

The only orbits which have their real turning points exactly at the focal plane 
are the periodic orbits. A periodic orbit as displayed in Fig. [To] (a) generates 
exactly two discrete points in the SOS, corresponding to the two distinct values 
of the radial distance p at which the axis z — is crossed. Both points in the SOS 
lie at p p = for the periodic orbit, as can be verified from the trajectory in the 
z - p plane which always crosses the z axis perpendicularly. The quasiperiodic 
trajectory of Fig. [l(] (b) corresponds to the leftmost curve in Fig. |ll|, which has 
its turning point almost at po in the SOS. The caustic is almost on the focal 
plane but still offset from it by an amount that is not discernible in Fig. [h] (b) . 

The distinction between the periodic orbit and its closely neighboring quasiperi- 
odic relative in the SOS of Fig. y] is appreciable - a pair of points generated by 
the former, versus a one-dimensional curve for the latter. But exactly on the 
line p p = 0, there exists an infinite number of other pairs of points, belonging 
to the periodic orbits of the same family. 



34 



_1 o ■ aa 0.4 o.e OS 1 



Figure 11: Poincare surface of section of the ray dynamics for L z = 0.1 fas in 
Fig. ^ (a-c). The minimum distance po from the z axis, given by Eq. (87), is 
indicated by the dashed line. All quasiperiodic trajectories fill smooth curves 
extending between some minimum p m i n > po and p = 2f. The two gray dots 
are the crossings of the focal plane generated by a periodic orbit as shown in 
Fig. M (a). The plot uses 2 / as the length unit. 



The SOS in these coordinates allows us to see directly in which places the 
focal plane comes into contact with the rays under consideration. This is a 
central piece of information when it comes to estimating the focussing at this 
plane where the quantum well is assumed to be. The forbidden regions around 
the z axis induced by the angular momentum barrier show up as inaccessible 
portions of the SOS toward small p. 

6.6 Accuracy of the semiclassical approximation 

We can also comment on the striking accuracy of the semiclassical approach 
in this system. The Poincare section shows that almost all trajectories (with 
the exception of the periodic paths) generate curves with the same topology: 
they begin and end at p = 2 / , with one turning point inbetween. There are 
other integrable systems for which the Poincare section has a more complicated 
structure, one closely related example being the ellipsoidal cavity or its two- 
dimensional counterpart, the ellipse billiard |2j|. In that case, the phase space 
consists of two components in which the topology of the trajectories is different: 
One type of motion consists of rays circulating around the perimeter as so- 
called whispering-gallery orbits, the other is a bouncing-ball oscillation across 
the short diameter J2(|. There is a division between these two types of trajecto- 
ries, similar to that between oscillation and rotation in a pendulum - called the 
separatrix. The WKB approximation or its higher-dimensional generalization, 
named after Einstein, Brillouin and Keller (EBK), cannot be applied without 
severe corrections in the vicinity of such a separatrix in phase space ||^, . 
In our case, this breakdown never occurs, and semiclassical results are thus of 
high accuracy. Being a conic section, the parabola can of course be considered 
as a limiting case of the ellipse, with one of its foci moved to infinity. This 
leaves no possibility for bouncing-ball trajectories, which leads to the absence 
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of a separatrix. 

Finally, it is worth asking why the sharp corners at the intersection between 
the paraboloid and the focal plane do not cause any corrections to our semi- 
classical treatment, even though the surface curvature at these points is clearly 
much shorter than the wavelength. It is known that in such cases diffraction can 
occur which makes it impossible to explain the mode structure purely based on 
classical orbits^, |24|. However, this phenomenon is absent for certain special 
angles subtended by the corners. One of these "benign" angles is precisely the 
90° angle we encounter at the corners of the double paraboloid, cf. Fig. ||. When 
the confocal condition is violated so that deviations from a right angle occur at 
the corners, we have to expect diffractive corrections to the semiclassical anal- 
ysis, resulting from classical rays that hit the corners and are reflected in an 
arbitrary direction because the law of specular reflection is undefined in that 
instance. Fortunately, we shall see in Section fj] that such orbits are far removed 
from the regions of phase space where we expect the important focussing modes 
to lie. 

In this section, we have discussed how the ray dynamics develops caustic 
structure, and how the latter can be represented with the help of the Poincare 
section. We have also observed that the high-intensity regions in the wave 
solutions correspond to the ray caustics, because there the density of rays is 
high - in fact divergent if we recall the discussion of the classical turning points 
in the effective potential below Eq. (E^). Therefore, even in situations where we 
cannot obtain the wave solutions easily, their possible intensity distribution can 
be inferred by investigating the ray dynamics first. This will now be carried out 
for a cavity that deviates from the ideal model shape. 

7 The non-confocal double paraboloid 

Having obtained an overview of the types of ray motion that can be encounterd 
in the parabolic dome, and established the connection to the mode structure 
of the full vectorial wave equations via the short-wavelength approximation, we 
now want to introduce a model cavity for which the wave solutions cannot be 
obtained by separation of variables. The variety of possible deviations from the 
ideal model geometry of Fig. || is enormous, so we have to restrict attention to 
certain special distortions that can be expected to be generic in some sense. 

7.1 The model deformation 

The distortions we choose are obtained by pulling the two intersecting paraboloids 
in Fig. [3] apart or pushing them together along the z axis by an amount 2 e. 
Specifically, in spherical coordinates as a function of polar angle 6, the shape is 
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given by 

r(g) = 2/ „ + ; 26 =. (91) 

l + costf 1 + ^1 + 6(1 -cos 2 61) 

The respective foci, which coincide in the integrable model, then move off the 
p axis. This non-confocal arrangement of the parabolic walls can be viewed as 
a model for fabrication-induced deviations from the ideal cavity shape - where 
the dome could be slightly too flat or too pointed. It can also be interpreted in 
a different way, taking into account the possibility that the boundary condition 
at the base of the dome is not exactly given by Eq. (Q), if some penetration 
of the field through the dielectric mirror on the quantum well is taken into 
account. This is of course a realistic expectation, and its effect on the wave 
solutions would be that the TE electric field no longer needs to be strictly 
symmetric under reflection at the focal plane. If one maintains that the dome 
has indeed been fabricated with its base in the focal plane, this "soft" boundary 
condition on the mirror can be modeled by assuming that our solutions should 
correspond to waves reflected at a plane removed from the dielectric interface 
by some amount e. 

Therefore, the non-confocal double paraboloid is a way of taking into account 
the cumulative effects of fabrication uncertainty and soft boundary conditions at 
the dielectric mirror with a single model parameter e, denoting half the distance 
between the foci of the top and bottom parabolic wall in the unfolded cavity. 
One could think that a perturbation theory in e could allow us to use the solu- 
tions obtained so far and smoothly extend them to the non-confocal situation. 
This is the traditional approach in physics and it is the reason why only simple, 
integrable systems are treated in textbooks on quantum mechanics or classi- 
cal mechanics alike. However, perturbation approaches become tedious and 
even impossible for wave equations whose short-wavelength limit (i.e. ray pic- 
ture) exhibits chaotic dynamics. The difficulties that arise can already be seen 
without introducing chaos, if we try to obtain the wave functions of an ellpsoid- 
shaped resonator as a perturbative expansion starting from the eigenfunctions 
of a spherical cavity. This poses no problems as long as one is interested only 
in modes of the ellipsoid whose topology is analogous to that found in the circle 
[ p8| . However, as mentioned earlier, the ellipsoid exhibits separatrix structure 
in phase space because there exists a type of motion that the sphere does not 
possess: the bouncing-ball trajectories. 

Analogous nonperturbative effects arise in the present model, because the 
distortion can lead to new types of trajectories that are not present in the con- 



focal cavity, in a process known as bifurcation [129, 30 . The first consequence of 



the deformation e is that the infinite families of periodic orbits break up, leaving 
only a distrete number of periodic orbits of the same topology, which can be 
divided in an equal number of stable and unstable paths. Stable paths have the 
property that rays with slightly different initial conditions remain close to the 
given periodic path for all times, while unstable periodic orbits are surrounded 
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Figure 12: Surface of section at L z = 0.1 / of a non-confocal double-paraboloid, 
with foci pulled apart by e = 0.02/. This destabilizes the cavity, leading to 
chaotic ray dynamics which generates an irregular cloud of points filling almost 
the whole region that is accessible for this L z . A special point is encountered 
on the line p p = where the irregularity gives way to a confluence of hyperbolic 
traces whose vertices are centered on a single, unstable periodic orbit. The 
spatial pattern of this new periodic orbit is shown in Fig. 13 (a). 



in their immediate neighborhood by chaos - trajectories deviate from such a 
periodic orbit at an exponential rate if the initial condition is only infinitesi- 
mally varied. For more quantitative statements and further background on the 
transition to chaos, the reader is referred to the literature [|9[ [[l| \n\ . 



7.2 Unstable and stable ray motion in the deformed cavity 

The Poincare section is very suitable as a diagnostic tool to identify this process 
of emerging chaos on on hand, and the stabilization of certain periodic orbits 
on the other hand. This is illustrated in Fig. [l2|. The perturbation consists of 
pulling the foci of the walls apart by e = 0.02 / along the z - axis. Since this 
preserves the axial symmetry of the cavity, L z is still a conserved quantity - 
the arguments of Section |6.3| rely on no other symmetries of the problem. We 
chose L z = 0.1/ in the plot. The small distortion of one percent is already 
sufficient to change the phase space portrait significantly, compared to Fig. 
Til). The unstable periodic orbit appearing prominently in Fig. 1^) as a so- 



called hyperbolic point, is shown in its spatial pattern in Fig. |13| (a). It is a 
self-retracing periodic orbit because it reflects from the boundary at normal 
incidence (in the z - p plane). 

The effects that chaos can have on the ray motion are illustrated in Fig. 
[l3| (b). Shown there is a single ray trajectory which superficially has some 
similarity to Fig. |l^ (b). However, the path does not trace out a well-defined 
caustic in Fig. O (b). What looks like a caustic here is in fact better described 



as two caustics of the type in Fig. 10, arranged almost symmetrically with 



respect to the focal plane. Note in particular the symmetric occurrence of cusps 
both below and above the line z = 0. Recall that in the integrable case the 
position po > z o of the caustic singularity is uniquely given by the turning points 
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Figure 13: Trajectories in the non-confocal cavity, (a) shows the unstable peri- 
odic orbit arising at e = 0.02 /, L z — 0.1 /, cf. Fig. [12| For the same parameter, 
a chaotic trajectory is seen in (b). Oscillatory motion around stable periodic 
orbits occurs in (c) and (d), where e = —0.02/ and L z = 0.1/ as in the SOS 
of Fig. O. The patterns of type (a) and (c) derive from the periodic motion 
of Fig. [nj (a) as a result of the shape perturbation. In real three-dimensional 
cartesian space, (e) shows the straight-line ray motion (arrows) giving rise to 
the curved "envelope" surface whose cross section we see in (a). 



£o and rjo (or equivalently uq, Vq), in the effective potential. Reversing the sign 
of the cusp coordinate zq corresponds to exchanging the role of £ and 77. The 
significance of Fig. [l^ (b) is therefore that the degrees of freedom £ and r\ are 
no longer decoupled, because during a single ray trajectory both the cusps at zq 
and — z are reached. By virtue of Eq. (|9|), the quantity Z is thus not conserved 
anymore. A trajectory is able to exhibit multiple points of closest approach to 
the z axis and is not strictly guided by caustics. 

Under these circumstances, it is not clear what to expect for the mode struc- 
ture of the cavity because we lose the possibility of labeling the eigenstates by 
a complete set of quantum numbers. This does not imply there arc no modes 
associated with chaotic rays, but one requires additional techniques to perform 
a semiclassical quantization |n , 32 . The destruction of the conserved quantity 
Z means that there is one less constraint which the ray trajectories have to sat- 
isfy; this allows them to fill two-dimensional areas instead of one-dimensional 
curves in the SOS. Since the SOS gives us a picture of how the rays intersect 
the plane z = 0, chaotic rays can be seen to show less concentrated overlap with 
that plane. We anticipate that the presence of true caustics is required to create 
the best focussing action. With this hypothesis, the goal must be to identify 
ray orbits that exhibit caustics. This occurs in the vicinity of stable periodic 
orbits, due to the fact that perturbed trajectories execute an oscillatory and in 
general quasi-periodic motion around such stable orbits. In Fig. |l2|, however, 
no stable periodic orbits can be identified, telling us that for the deformation 
chosen ther, no stable modes with L z = 0.1 / should exist. 

The situation changes if we consider Fig. [l4|, in which L z is the same but 
the sign of the non-confocal displacement e is reversed. The walls of the dou- 
ble paraboloid are hence pushed together instead of being pulled apart. The 
resulting phase space structure in the SOS differs markedly from Fig. ^2[ many 
trajectories trace out one-dimensional curves in the SOS which organize as closed 
loops, forming island chains that proliferate with various sizes. All these islands 
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Figure 14: Surface of section at L z = 0.1 / with e = —0.02 /. The cavity devel- 
ops a multitude of stable periodic orbits surrounded by elliptical islands. The 
most prominent island of stability is centerd on the the line p p = where a small 
circle indicates the location of the corresponding periodic orbit around which 
other trajectories can oscillate. Shown in the inset is the central stable peri- 
odic orbit. The next innermost closed line in the SOS belongs to the trajectory 
shown in Fig. |l3| (c). 

are centered around stable periodic orbits - the biggest island of stability lies 
symmetrically around the line z = and corresponds to oscillatory motion of 
the type shown in Fig. [l3| (c). The center of the island is in fact formed by 
a periodic orbit similar to Fig. [l3| (a) - the only difference being, that small 
perturbations of its initial conditions do not lead to chaos as in Fig. [lj, but to 
the motion of Fig. |l^ (c). 

Another oscillatory ray path centered at a stable periodic orbit is shown in 
Fig. |l3| (d). The pattern should be compared to Fig. |l0| (c) which has the same 
L z . The similarity is apparent, except for the fact that the path in Fig. |l^ (c) 
will eventually fill the remaining gaps in that plot, if one follows it longer. The 
path in Fig. |l3| (d), on the other hand, is truly restricted to the vicinity of a 
self-retracing orbit which reverses its propagation direction at one end due to 
perpendicular reflection at the wall, and at the other end by running up the 
centrifugal barrier perpendicular to the z axis. 

All islands of stability in Fig. [l4| generate their own caustics, which are topo- 
logically different from the ones in the integrable system. The caustic created by 
the orbit in Fig. [l^ (c) is simply the boundary of the regions into which the ray 
never penetrates. The difference between the absence and presence of caustics 
in Figs. |l^ (b) and (c) is not easily appreciated if we consider only the real-space 
plots. Here, the usefulness of the Poincare section as a diagnostic tool is again 
to be noted - showing two-dimensional clouds of points versus one-dimensional 
curves, respectively, for trajectories without and with caustics. 

It follows from the preceding discussion that a negative e leaves us with a 
cavity that is in many respects similar to the unperturbed double paraboloid, cf. 
the ray pattern of Fig. [l^ (d). However, the qualitative and important difference 
is that some periodic orbits are now more stabilized than at e = 0. In particular, 
there are simple ray bundles such as Fig. [l3] (c) that promise reasonable focusing 
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Figure 15: For the same deformation e = 0.02 / as in Fig. [lj, this surface of 
section at the smaller angular momentum L z = 0.03 / shows that a stable orbit 
exists in addition to the unstable hyperbolic one. This is indicated by the elliptic 
(lens-shaped) island structure. The hyperbolic point is located to the right of 
the island. The ray pattern near the stable periodic orbit is shown in the inset. 
The small corresponding mode volume is apparent. 



close to the center of the unfolded cavity. The physical explanation for the 
general stabilizing effect that we achieved by moving the paraboloids closer 
together lies in the well-known fact that a two-mirror resonator configuration 
has a a focusing action when the mirrors are separated less than the sum of 
their radii of curvature. Conversely, mirrors that are further apart than this 
criterium act in a defocusing way. This is consistent with the observation of a 
large chaotic domain in Fig. [l2| 

These simple arguments, and the chaotic picture of Fig. |l2|, seem to suggest 
that stable ray motion is not to be expected in the supposedly defocusing con- 
figuration with e = 0.02/. However, when applying the standard criteria for 
focusing and defocusing resonator geometries, we have to bear in mind that we 
are dealing with a centrifugal billiard whose ray trajectories are curved. The 
effect of the centrifugal barrier is to push the regions of allowed ray motion out- 
wards until only a small patch surrounding the equatorial corners of the cavity 
is accessible. At large L z the motion is then so confined that chaos does not de- 
velop. This is just the whispering-gallery phenomenon JT^] . On the other hand, 
at L z = 0.1 / we certainly found chaos with no remaining islands of stability. 
Small L z are what we must be interested in if concentration near the focal points 
is to be achieved. 

In view of this, it is all the more surprising that the same cavity does in fact 
support stable orbits at even smaller angular momenta than in Fig. [l2|. This is 
shown in Fig. \LB for L z = 0.03 /. The periodic orbit responsible for the single 



stable island in that SOS is again almost identical to the one shown in Fig. 13 
(a), and its oscillatory neighborhood is analogous to Fig. |l3| (c); the inset of 
Fig. [ll] shows this similarity. This stable orbit exists only at sufficiently small 
L z \ its associated island in the SOS shrinks to a point when L z w 0.038 /. The 
conclusion is that both the nominally focusing and defocusing configurations 
e = ±0.02 / permit the formation of ray bundles with a spatial distribution as 
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in Fig. [T^ (c), and hence the stable modes associated with this pattern should 
be robust. This is also confirmed by analogous Poincare sections for larger 
displacements of the foci. At larger |e|, the motion of type Fig. |l3| (c) and the 
inset of Fig. [15] is in fact stabilized further - for both directions of displacements 
alike. 

The modes corresponding to this particular ray pattern are closely related 



to the fundamental s-waves we discussed in Section 3.5, because both arise from 
ray bundles in the immediate vicinity of the shortest periodic orbits in the cavity. 
For e = this was the family of paths in Fig. [To] (a), members of which can be 
smoothly deformed into Fig. [l3| (a) without changing the topology - i.e., the 
number and sequence of reflections and turning points. We shall therefore call 
all these orbits the fundamental orbits of the cavity. The mode spacing of the 
corresponding eigenstates should be comparable as well for the perturbed and 
unperturbed case. However, we have to defer a detailed analysis of the wave 
solutions and their semiclassical correspondence to a future paper. Here, the 
goal has been to introduce the ray dynamics and its phase space as the backbone 
on which the mode structure is built. 

Assuming that the deformation is e = 0.02/, we have the peculiar situation 



that the fundamental orbit is unstable if L z > 0.038/, cf. Fig. 12. Therefore, 
the most desirable modes will be those with smaller L z . According to Eq. 
we have to choose modes with low n and high k to achieve this. For the 



experimental cavity we have kf w 14tt. Taking n — 1 as in Section 3.5, we 
arrive at the semiclassical value 

L* = tj-« 0.023/ (92) 

14 7T 

which is close to the situation depicted in Fig. O. The difference in the SOS 
is insignificant. We have no accurate way of determining the actual value of 
e most closely describing the real structure, but these considerations give us 
considerable confidence that modes with a spatial pattern as in Fig. |l3| (c) or 
Fig. [l^ will be found in the cavity, because the relevant L z estimated above is 
in a range where this fundamental orbit is stable - irrespective of the sign of e 
and moreover largely independent of its magnitude. 



8 Bragg mirror as an escape window in phase 
space 

The internal ray dynamics of the dome resonator has up to this point been 
evaluated under the assumption that the cavity is a perfect resonator. There 
are two physical mechanisms that invalidate this viewpoint: absorption in the 
gold mirror and transmission through the Bragg grating. The trade-off between 
the comparatively large absorption of a metal on the one hand and its ability 
to reflect omnidirectionally have been discussed in Ref. [B3|. In our context, 
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metallic absorption will always degrade the Q factor because the gold layer pro- 
vides only an estimated 95% reflectivity J33|. However, the reflectivity of the 
Bragg mirror can be significantly lower for certain modes and in that case con- 
stitutes the dominant mechanism for Q-spoiling. The variable that determines 
the reflectivity of the Bragg mirror (at the fixed operating frequency) is the an- 
gle of incidence x with respect to the z-axis. For the purposes of a qualitative 
analysis, we assume that the Bragg reflectivity is unity for x < 22° = Xc but 
drops to sw 20% outside this cone of incidence In other words, Xc is the 

boundary between absorption-dominated and leakage-dominated Q factors. A 
second window of high reflectivity opens for rays at very oblique incidence on 
the grating surface, more specifically for x > 60°. This second window will be 
discussed further below. 

The ray picture allows us to use this rough transmission criterion as a guide 
in order to separate long-lived cavity modes from short-lived ones. The angle x 
between z-axis and a trajectory is, according to Eq. (|72|), given by 

cosx = r ■ e z = z, (93) 

so that 

f - z 2 = sin 2 x- (94) 



One can substitute this as the righthand side of Eq. (|82|) and obtains an equation 
for a curve in the plane p, p (p = p p ) spanning the Poincare section: 



M = ysin 2 X -^f. (95) 

Using the critical value of Xc in this equation specifies the escape condition in 
the Poincare section: the Bragg mirror becomes ineffective when 



bp I > ysin 2 Xc -^|. (96) 

In order to get a feeling for the type of ray orbits that can remain in the 
cavity under this escape condition, we plot in Fig. ^| the resulting curves in 
the surface of section for the two different values of L z appearing in Figs, [n], 
p"2| , pT| and |f5|. The plot should be superimposed on these plots to decide which 
parts of the respective phase space falls within the high-reflectivity range of the 
DBR grating. Note that the critical lines for ray escape are independent of 
deformation because they rely only on Eq. (|S3|). 

As a result of this comparison, we find first of all that low angular momenta 
are required by the escape criterion, because the phase-space area enclosed by 
the critical curves in Fig. shrinks with increasing L z . This is understand- 
able because the ray motion in this case has a strong azimuthal component 
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Figure 16: Assuming that the Bragg mirror stack at the plane z = yields 
high reflectivity only for waves within \c = 22° from the surface normal or for, 
the regions of high reflectivity in the Poincare section are bounded by Eq. ( p6| ) . 
The resulting curves bounded by \p p < 0.374 are shown for L z = 0.03 (solid 
line) and L z = 0.1 (dashed). A second high-reflectivity window exists for rays 
falling between the boundary of this plot and the solid line near the boundary. It 
becomes relevant only for the integrable confocal cavity because the perturbed 
shapes have no stable orbits in this second window. 



contributing to the tilt angle with respect to the z-axis. Let us turn our atten- 
tion to the stable periodic orbits arising in the chaotic Poincare sections. The 
case L z = 0.1 / shown previously for illustrative purposes turns out now to be 



roughly the maximum angular momentum at which the stable orbit of Fig. 14 
is still confined by Bragg reflection. The lower angula momentum L z = 0.03 / 
coming close to the estimated value for the s- waves of our experimental cavity, on 
the other hand, places the stable periodic orbit well inside the high-reflectivity 
range of the DBR. For the case of a defocusing deformation this is illustrated 
in Fig. H[ The peridic point is at p w 0.086. For a focusing deformation of the 
same magnitude, e = —0.02, the periodic point lies at p sa 0.99. Both values are 
to the right of the solid line in Fig. [l6], corresponding to high reflectivity. 

For the chaotic orbits, we observe that they spread out over the Poincare 
section in such a way as to yield significant overlap with the low-reflectivity 



regions of Fig. 16, This is true for all Poincare sections shown in this paper. 
Therefore, we conclude that cavity modes associated with the chaotic phase space 
regions are short-lived, and the corresponding broad resonances will not affect 
the spontaneous emission enhancement of the parabolic dome. A quantitative 
estimate of the resonance lifetimes could be obtained by measuring the time that 
a chaotic trajectory spends, on average, in the high-reflectivity region without 
excursions beyond the critical line. However, we shall not attempt quantitative 
predictions at this stage of our investigation, and defer it to future work. 

A quantitative analysis would also be necessary to determine the modal life- 
times in the marginal case of the ideal confocal cavity. The reason is that the 
ray picture alone does not allow a clear distinction between classically confined 
and unconfined orbits, because the classification according to stable and unsta- 
ble trajectories does not apply in the integrable parabolic dome. All the solid 
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Figure 17: The Poincare section combined with the escape conditions can be 
used to extract information about the lifetime and escape locations. This is 
illustrated here for a single ray orbit (black trace) , followed for 500 crossings of 
the focal plane. The gray area is the region which has to be avoided by the ray 
in order to remain in the cavity. 



curves in the Poincare section of Fig. [Tl] cross into the low-reflectivity region 
of Fig. [If] at some point, but the time spent in the high-reflectivity range can 
be very long classically. To illustrate this, we show in Fig. [Tt] a particular ray 
trajectory for L z = 0.1 / in the confocal paraboliod, which for almost 500 cross- 
ings of the focal plane remains inside the regions of high reflectivity. This time, 
the second window of high reflectivity close to the border of the SOS is im- 
portant because the ray alternates between the low- and high-x windows from 
one crossing of the focal plane to the next. The regular nature of this motion 
makes long lifetimes possible because it strictly prevents the ray from entering 
the low-reflectivity region for long times, whereas a chaotic orbit would quickly 
explore this domain in a quasi-random way. 

The trajectory shown in Fig. [l?] is practically identical to the one shown in 
Fig. nil (b). The alternating way of intersecting the focal plane can be under- 



stood from that figure, or from Fig. 10 (a) which shows periodic orbits closely 



neighboring the quasiperiodic trajectory of plot (b). Note that the ray model al- 
lows us in addition to predict the spatial location where the mode corresponding 
to this ray bundle will preferentially be coupled out through the Bragg mirror. 
As can be seen in Fig. [l7], the low-reflectivity region is reached for the first 
time when, after many reflections, the trajectory departs from the immediate 
neighborhood of the focal region, i.e. intersects the focal plane with a p that is 
slightly too large. 

The subtle balance of parameters that prevents chaos from appearing will, in 
all experimental realizations, be shifted to either the defocusing or the focusing 
side. Therefore, the above ray analysis of the mixed phase spaces for these 
two situations above is our main concern. However, as in the previous sections 
the integrable case is a useful starting point to illustrate our strategy. The 
advantage of the ray approach is that it provides fast and intuitive predictions, 
but further studies are required in order to determine how this model succeeds in 
characterizing the cavity quantitatively. Paradoxically, we can already conclude 
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that the existence of chaos and islands of stability makes it easier to obtain 
results from a ray analysis, because there is a sharper separation between long 
lifetimes for the stable modes discussed above and short lifetimes for modes 
associated with the chaotic portions of the SOS. 

9 Conclusion 

In this paper we have examined the modal structure of the electromagnetic field 
in a semiconfocal piano-parabolic cavity (or, equivalently, in a double-paraboloid 
confocal cavity) in view of our recent fabrication of semiconductor microcavities 
having that geometry. In order to account for the effects of the inevitable 
fabrication defects we also considered the stability of the modes with respect to 
deformations consisting of deviations with respect confocality. This theoretical 
analysis was thus motivated by our ongoing experiments on these structures, 
and feeds back into this experimental work by opening a novel perspective in 
terms of investigating the chaotic structure and dynamics of some of the modes 
of cavity. 

Regarding the structure of the modes in the parabolic cavity, we note that 
the scalar wave equation is solvable analytically by separation of variables. How- 
ever, the vectorial boundary conditions for the electromagnetic field destroy this 
property, leaving only the cylindrical symmetry. Nevertheless, the fundamental 
series of s-waves (free of azimuthal nodes) in a confocal electromagnetic cav- 
ity can be solved rigorously. It has its energy concentrated in a small volume 
(of order A 3 ) around the focal point, even though at the focal point itself the 
electric field is zero due to the vectorial nature of the field. The higher order 
modes cannot be solved as readily in the full three-dimensional model, but it is 
possible to appreciate their features by reducing the problem to scalar form. In 
these higher order modes, the energy is concentrated in lobes that surround the 
focal point but avoid it because of the centrifugal barrier that arises from the 
cylindrical symmetry. Indeed, these modes correspond to non-zero values of the 
angular momentum (m ^ 0) and for large values of m tend towards a type of 
whispering-gallery modes with intensity concentrated in a ring along the focal 
plane, [cf. Fig. |o| (d)]. 

The stability of the modes of the parabolic cavity with respect to geometrical 
deformations can be assessed by examining the ray trajectories that correspond 
to each mode. For a deformation that corresponds to a small deviation from 
confocality, chaotic ray patterns emerge. However, we also find stable ray orbits 
concentrated in a small part of the cavity volume. Independent of deformation, 
the most important stable orbits being those which in cylinder coordinates p 
and z follow the shortest possible periodic trajectory. This general topology is 
the same for a range of deformations (including the ideal confocal cavity) and 
corresponds to a ray returning to the same p and z after two reflections, missing 
the focal point by a small amount because the field there has to vanish. The 
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generic shape of this orbit is represented in Fig.|l3], and its special modification in 
the confocal case with its marginal stability is shown in Fig.|l^. The topological 
equivalence between the stable orbits of the distorted cavity on one hand and of 
the confocal system on the other indicates that the structure of the fundamental 
s-wave is stable with respect to deformations. 

From the experimental viewpoint, the results of this theoretical analysis indi- 
cate that the cavities already fabricated in our laboratory should possess stable 
modes in which the energy is confined in a volume of order A 3 in the vicinity 
of the focal point, in spite of fabrication errors. The higher order modes, in 
which the field is concentrated away from the focal point, in whispering-gallery 
type configurations, will be unstable because of the presence of fabrication de- 
fects. At the same time these modes will decay very fast as they correspond to 
oblique incidences onto the Bragg mirror, at angles for which the mirror is no 
longer reflecting. Experiments are in progress to characterize the structure and 
dynamics of both the stable and unstable modes . The robust stable modes 
in which the field is confined in the vicinity of the focal point should give rise 
to strong enhancement of the spontaneous emission of a dipole (such as a semi- 
conductor quantum well or a semiconductor quantum box) placed there, and 
a concomitant lowering of the lasing threshold, even for our cavities that are 
of mesoscopic dimensions. This is because because even in such large cavities, 
whose geometric volume is of the order of a few thousand cubic wavelengths, 
the central lobe of the fundamental s-wave (which contains most of the energy) 
has an effective volume of the order of one cubic wavelength. 

These considerations underscore the interest that parabolic microresonators 
present by exhibiting quantum electrodynamic effects as well as optical chaos, 
in spite of their relatively large dimensions. In addition, the mesoscopic cav- 
ity dimensions of these structures are an important practical feature, as they 
make the fabrication accessible to existing experimental techniques (such as Fo- 
cused Ion Beam etching) while, at the same time, they greatly facilitate the 
theoretical analysis of these devices as they permit the use of short-wavelength 
approximations . 
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A Parabolic coordinates 

The parabolic coordinates £, rj, <fi are related to the three-dimensional cartesian 
coordinates according to 




(97) 
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Or, equivalently, 




(98) 



where r = yfx 2 



y 2 + z 2 is the spherical radius vector. With this definition, £ 
and rj have the same dimensions as the cartesian coordinates, which is helpful 
for physical considerations. The surfaces £ = constant are paraboloids by revo- 
lution about the positive z-axis having their focal point at the origin, while the 
surfaces r\ = constant are directed along the negative £-axis. The plane z = 
corresponds to the condition £ = rj. In terms of the cylindrical coordinates 
■y 2 ,z,4> the parabolic coordinates obey 



V^ 2 



rj) 



(99) 



and 




(-VZ-i+y/V-v) 



(100) 



In these parabolic cordinates, the electric field E = (E$, E v , E<jj) is related to 
its representation in cylindrical coordinates according to 



E 




(101) 



I Et = ±(E + + E_) 
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